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ABSTRACT 

Context. TWA22 was initially regarded as a member of the TW Hydrae association (TWA). In addition to being one of the youngest 
(a 8 Myr) and nearest (» 20 pc) stars to Earth, TWA22 has proven to be very interesting after being resolved as a tight, very low- 
mass binary. This binary can serve as a very useful dynamical calibrator for pre-main sequence evolutionary models. However, its 
membership in the TWA has been recently questioned despite due to the lack of accurate kinematic measurements. 
Aims. Based on proper motion, radial velocity, and trigonometric parallax measurements, we aim here to re-analyze the membership 
of TWA22 to young, nearby associations. 

Methods. Using the ESO NTT/SUSI2 telescope, we observed TWA22 AB during 5 different observing runs over 1 .2 years to measure 
its trigonometric parallax and proper motion. This is a part of a larger project measuring trigonometric parallaxes and proper motions 
of most known TWA members at a sub-milliarcsec level. HARPS at the ESO 3.6m telescope was also used to measure the system's 
radial velocity over 2 years. 

Results. We report an absolute trigonometric parallax of TWA22 AB, n = 57.0 ± 0.7 mas, corresponding to a distance 17.5 ± 0.2 pc 
from Earth. Measured proper motions of TWA 22AB are ^„co.s(5) = -175.8 ± 0.8 mas/yr and yu^ = -21.3 ± 0.8 mas/yr. Finally, from 
HARPS measurements, we obtain a radial velocity Vraj = 14.8 ±2.1 km/s. 

Conclusions. A kinematic analysis of TWA22 AB space motion and position implies that a membership of TWA22 AB to known 
young, nearby associations can be excluded except for the y6 Pictoris and TW Hydrae associations. Membership probabilities based 
on the system's Galactic space motion and/or the trace-back technique support a higher chance of being a member to the /? Pictoris 
association. Membership of TWA22 in the TWA cannot be fully excluded because of large uncertainties in parallax measurements 
and radial velocities and to the uncertain internal velocity dispersion of its members. 



Key words. Astrometry: trigonometric parallaxes - 
~ Galaxy: open cluster and association 



Techniques: radial velocities ~ Stars: binaries, distances, fundamental parameters 
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1. Introduction 

Over the past decade, various coeval moving groups of stars 
have been discovered in the solar neighborhood, includ- 
ing three well-known associations, TW Hydrae (hereafter 
TWA, [Ka stner et al. 1997 ), /3 Pictoris (hereafter /3 Pic, 
[Barrado yNavascues et al. 1999| l, and Tucana/Horologium 
(hereafter Tuc-Hor, IZuckerman & Webb 20001 

[Torres et al. 20001 1. These groups occupy an important as- 
trophysical niche thanks to their proximity (< 100 pc) and youth 
(< 100 Myr). They offer the best targets for various studies, such 
as imaging searches for young brown dwarf and planetary mass 
companions, proto-planetary or debris disk-related programs, 
etc. Any young (sub)stellar binary members in tight orbits 
in these moving groups can serve as valuable calibrators for 
evolutionary models. 



to conclude that TWA22 could indeed be the nearest TWA 



Song et al. (2003) identified the TWA22 system as a new 



member of TWA based essentially on the presence of a strong 
Li /16708 absorption feature and its close proximity to TWA in 
the sky. A further proper motion analysis led |Scholz et al. (2005)| 

Send offprint requests to: teixeira@astro.iag.usp.br 

* Based on observations performed at the European Southern 
Observatory, Chile (76.C-0543, 077.C-01 12, 078.C-0158, 079.C-0229). 
** Table 4 is only available in the electronic edition of the journal. 



member to Earth. However, more recently, Mamajek (2005) 
has performed a convergent point analysis using several TWA 
members. He found that TWA22 has a low probability of 
membership in the TWA. This conclusion was refuted by 



Song et al. (2006) arguing a lack of reliable distance determi- 



nations for most TWA members to firmly reject TWA22's mem- 
bership. ISongjetair(2006)jargued instead that the strong lithium 
line seen at TWA22, rarely seen among other stars with similar 
spectral types, implies a probable membership to TWA or /3 Pic. 

During a VLT/NACO deep-imaging survey for close com- 
panions to stars in young associations, TWA22 was resolved as 
a tight (a: 100 mas) binary with a projected physical separation 
of 1.76 ± 0.10 AU (see Bonnefoy et al. 2009). 

Regardless of TWA22's membership in TWA or /3 Pic, 
TWA22 AB is a precious dynamical mass calibrator because of 
its young age (a: 10 Myr) and the low masses of its components. 
A well-known age and a distance are crucial to the validity of 
calibrating evolutionary model calculations. Therefore, there is 
need for accurate observational data, such as the trigonometric 
parallax, proper motions, and radial velocity. 

As a part of a larger project for the trigonometric parallax de- 
termination of all TWA members, we performed astrometric and 
photometric observations of TWA22AB at ESO NTT/SUSI2 
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(La Silla - Chile). Radial velocity measurements were also ob- 
tained to complete the kinematic analysis of TWA22 AB. Using 
these new data, we performed a kinematic analysis to test the 
membership of TWA22 AB in TWA, jS Pic, and Tuc-Hor. Our 
analysis uses only moving- group members with Hipparcos mea- 
sured distances. In addition to Hipparcos moving -group mem- 
bers, we used 2M1207A whose accurate trigonometric distance 
was obtained from gound-based measurement (Ducourant et al. 
2008). The possibility of TWA22 AB being a member of other 
nearby groups was investigated and discarded by more evident 
incompatibilities. 

We present observation and data reduction in Sect. 2, and in 
Sect. 3, we discuss the membership of TWA22 in TWA, [3 Pic, 
and Tuc-Hor The conclusion is provided in Sect. 4. 

2. Observations and data treatment 

2.1. Trigonometric parallax and proper motion 

We used ESO NTT-SUSI2 which provides a good compromise 
between a large field of view (5.5' X 5.5') for a sufficient sam- 
pling of background stars and a small pixel size (80.5 mas) nec- 
essary for sub-milli arcsecond astrometry. Data were acquired 
over five observational epochs, and all observations were ob- 
tained around the meridian transit (hour angle < 0.5h). In our 
imaging, we used an I-band filter to minimize the differential 
color refraction effects (DCR). Residual DCR effects were re- 
moved from single observations following the method described 
in [Ducourant et al. (2007)| Multiple images were obtained at 
each epoch to reduce astrometric errors. The image of highest 
quality is selected as a master frame that will be used for cross- 
identification, alignement and scaling. The alignment of CCD 
axes and the plate scale were estimated using 2MASS catalog 
sources (iCutri et al. 20031 1. 

A stellar point-spread function for each frame was fitted 
using the DAOPHOT II package (IStetson I987I I. Then, we cre- 
ated catalogs of measured positions (x, y), internal magnitudes, 
and associated errors for all stars on each frame. Observational 
data were processed through a global treatment, as described in 
Ducourant et al. (2007, 2008), and a solution was derived for 
TWA22 AB relative to background stars (14.5 < I < 18.5 mag). 
In our astrometry data reduction and analysis, we ignore any in- 
fluence of binarity. But, we assessed the effect of binarity in Sect. 

2.2. Then, a statistical conversion from relative to absolute par- 
allax and proper motions, based on the Besanfon Galaxy model 
jRobin etal. 2 003, 2004), was derived (An = 0.35 ± 0.01 mas, 
AlJ„cos(d) = -5.43 ± 0.02 mas/yr, Afig = +2.98 + 0.02 mas/yr). 
Both the final estimated TWA22AB [a = 10''17'«26.79', 6 = 
-53°54'26.5", epoch - 2006.763 yr) proper motion and trigono- 
metric parallax are given in Table [T] In a separate table (Table 
|2]), apparent and absolute Bes sel (V, R, I) and 2MASS (J, H, K) 
magnitudes ( ICutri et al. 2003l l are listed. 

Figure [T| presents the apparent displacement of TWA22 AB 
relative to the background stars due to the parallax and proper 
motion. 

Table 1. Absolute astrometric parameters for TWA22AB and 
radial velocity derived in this work. 



n 


d 






Vrad 


mas 


pc 


mas/yr 


mas/yr 


km/s 


57.0±0.7 


17.5±0.2 


-175.8±0.8 


-21.3±0.8 


14.8±2.1 



-100 -300 

Dx (mas) 

Fig. 1. Apparent astrometric displacement of TWA22AB to- 
gether with the best parallax and proper motion fit of the data. 
Blue dots correspond to weighted mean positions at each obser- 
vational epoch. 

Table 2. Apparent and absolute magnitudes of TWA22 AB: V, 
R, I Bessell filter data from this work and J, H, and data from 
2MASS. 



M 



(mag) 



V 


13.99 


±0.02 


12.77 


±0.19 


R 


12.50 


±0.14 


11.28 


±0.19 


I 


11.13 


±0.23 


9.91 


± 0.18 


J 


8.55 


±0.01 


7.33 


± 0.07 


H 


8.09 


± 0.04 


6.87 


± 0.10 


K 


7.69 


± 0.02 


6.47 


± 0.07 



2.2. Impact of binarity onto the parallax determination. 

With ESO NTT/SUSI2, it is not possible to resolve two com- 
ponents of TWA22 AB. Therefore, with our astrometric obser- 
vations we are measuring the photocenter of the system. This 
photocenter may not coincide with the center of mass because 
these two different positions depend on the ratio of mass and lu- 
minosity of the components. A photocenter of a binary should be 
affected by an elliptic, periodic movement with the same period 
as the binary orbit. The amplitude of this variation depends on 
the orbital parameters, mass, and luminosity ratios. 



Following I Van de Kamp (1967) our condition equations for 
TWA22AB were modified to include the photocentric move- 
ment. The condition equations are written for each star on each 
of the N frames considered (including the master frame). These 
equations relate the measured coordinates to the steUar astromet- 
ric parameters: 

Xo + AXq + jj.x(t - to) + nFxit) + aQa = aix(f) + a2y(t) -H 03 (1) 

Yo + AFo + Mt - to) + TTFrit) + aQg = bix(t) + b2y(t) + h (2) 

where (X^, Yq) are the known standard coordinate of the star at 
the central epoch fo and (x(f), y{t)) are the measured coordinates 
on the frame (epoch f) that need to be transformed into the mas- 
ter frame system. AXq, AYq, fix, I-Iy, ^, and a are the unknown 
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Table 3. Heliocentric coordinates (X, Y, Z) and Galactic velocities (U, V, W) of TWA22 AB together with mean values for TWA, fi 
Pic and Tuc-Hor 



Name 


X 


Y 


Z 


U 


V 


W 


Age 


d 






pc 






km/s 




Myr 


pc 


TWA22 


3.5 ±0.2 


-17.2±0.2 


0.7± 0.2 


-8.0±0.4 


-17.1±2.1 


-9.0+0.2 


<10(1) 


17.5+0.2 


TWA 


14.4±8.7 


-46.9±5.3 


22.5±2.5 


-10.0±2.0 


-17.6±1.4 


-4.8+1.1 


8(2) 


55+7 


jSPic 


9.2±28.9 


-7.5±13.0 


-13.2±6.5 


-10.9±1.9 


-16.1±1.0 


-9.0+1.2 


12(3) 


34+14 


Tuc-Hor 


14.1±17.7 


-19.4±8.0 


-34.9+3.5 


-9.5±1.7 


-20.6+1.7 


-0.6+2.6 


27 (4) 


46+5 



stellar astrometric parameters. Both AXq and AFq yield correc- 
tions of the standard coordinates of the star on the master frame, 
lix and are projected proper motions in right ascension and 
declination, n is the parallax, and a is a semi-major axis of 
the photocentric trajectory relative to a barycenter. Coefficients 
ifli, bi) are unknown frame parameters that describe the trans- 
formation to the master frame system, (F„, Fs) are the paral- 
lax factors, and (Qa, Qg) are known orbital factors (based on 
|Bonnefoy et al. (2009) orbital data). The unknown coefficient 
or IS given by a = a(R - j5) where a is the semi-major axis 
of TWA22 B's orbit around TWA22 A, R is a fractional mass, 
R - P is the fractional distance of the primary to the 

photocenter y6 = i+io-°-4a»i where Am is the magnitude difference 
between A and B. 

Although it is impossible to include the parameter a formally 
as a variable in our equations, we could externally determine 
its value to correct Xo, Fq for the orbital motion of photocenter 
around barycenter. For this, we assumed the mass-luminosity 
relations as: M oc (L)"^'^ and used the extrapolated magni- 
tudes in the I-band adapted from |Bonnefoy et al. (2009) Am/ - 
0.46 - 0.87 mag to determine aQa and aQs and correct Xq, Tq 
from these quantities. The resulting astrometric parameters for 
TWA22 AB are then: n = 56.4±0.7 mas, n„cos{6) = -178.2±0.7 
mas/yr, and - -9A + 0.8 mas/yr 

We note that the parallax value obtained from this study is 
within lcr„ of the one given in Table [T] where it did not con- 
sider the effect of photocenter' s periodic movement. On the other 
hand, we observe a large variation in the proper motion in dec- 
lination. This suggests that our astrometric data may not cover a 
long enough time interval to properly account for the 5.144yr, 
orbital periodic signal from ,Bonnefoy et al. (2009) The influ- 
ence of this signal over 1.2 yr (duration of the observational 
program of TWA22 AB) is mainly equivalent to an offset of the 
proper motion in declination. We therefore conservatively stick 
to the values given in Table[T]as the best fit to our data. 



2.3. Radial velocity 



TWA22 AB was observed with HARPS ( [Mayor et al. 2003| l 
over 2 years. HARPS is a high-resolution (R=l 15,000) fiber- 
fed cross-dispersed echelle spectrograph functioning on the 
ESO/3.6m telescope. For TWA22 AB (V = 13.99), we used 15 
min exposures to obtain 60 spectra with SNR ranging from 6 
to 10. The standard HARPS reduction pipeline was used to de- 
rive radial velocities from the cross-correlation of spectra with 
a mask for an M2 star. The instrument is generally stable over 
one night (nightly instrumental drifts < 1 m.s ' ), and we 
performed a precise nightly wavelength calibration from ThAr 
spectra ( [Lovis & Pepe IQOT ). Out of these measurements, we 
estimated a heliocentric radial velocity Vr - 14.8 ±2.1 km.s"' 
for the system. This radial velocity is likely to be affected by 
the SBl status of TWA22 AB. We took the binary nature into 



account in our uncertainty value by simulating expected ampli- 
tudes from the binary orbital motion based on the parameters 
from Bonnefoy et al. (2009). 

3. Kinematic analysis 

Based on new trigonometric parallax, proper motion, and radial 
velocity, we re-examine below the membership of TWA22 AB to 
TWA, yS Pic, and Tuc-Hor. Membership in other young, nearby 
associations was rejected by more evident discrepancies in age, 
distances, etc. 



3.1. Space motion 

To test the membership of TWA22 AB, our first approach is to 
statistically compare its galactic space motion with the mean 
UVW values for TWA, /? Pic, and Tuc-Hor members. Only stars 
with known trigonometric parallaxes were considered here be- 
cause accurate distances are crucial in the UVW calculation. 
For TWA, our sample set includes TWA 01, TWA 04, TWA 09, 
and TWA 1 1 with trigonometric parallaxes from the Hipparcos 
catalog ( ESA 1997) and 2M1207A from |Ducourant et al. (200"8)1 
Although TWA 19 has a Hipparcos parallax, we did not consider 
it here because Mamajek (2005) classified it as non TWA mem- 
ber SSSPMJl 102-3431 TWA member (ITeixeira et al. 200 8 ) has 
also been excluded from our study because no radial velocity 
was available. Their radial velocities were extracted from Table 1 
of Mamajek (2005) (T orres et al. 2003 data for T WA 01 and 
TWA9A; Tor res etal. 19 95 for TWA4: lReid2003] for TWAll; 
and Mohanty etaL2003jfor 2 M1207A). For jS Pic we used the 
list of suggested members by [Torres et al. (2006)| Astrometric 
measurements and radial velocities were respectively obtained 
from the Hipparcos catalog dESA 1997] l and the Table 6 of 
[Torres et al. (20067[ Finally, for Tuc-Hor, astrometric measure- 
ments and radial velocities were respectively obtained from 
the Hipparcos catalog dESA 199711 and K harch enko et al. (2007)[ 
for suggested members by Zuckerman & Song (2004)[ and 
[Torres et al. (2008), The selection of stars and the original data 
used in our analysis are presented in Table |4] (only available 
in the electronic edition). Calculated mean space motions of 
TWA, j6 Pic, and Tuc-Hor are reported in Table |3] together 
with their spatial heliocentric coordinates and velocities for 
TWA22AB. In this Table positive X(U) points to the Galactic 
center, Y(V) is positive in the direction of Galactic rotation 
and Z(W) is positive toward the north Galactic pole. Ages 
are from (l)rSong et al. (2003)[ (2)|De la Reza et al. (2()06 )1 (3) 



[Ortega et al. (2004)[ (4) [Makarov (2007)[ " 

The mean values derived here (Table [3]l for TWA, /3 Pic, 
and Tuc-Hor are in good agreement with published values : (- 
11, -18, -5), (-11, -16, -9), and (-11, -21, 0) from Table 7 of 



Zuckerman & Song (2004)[ and (-10.5 ±0.9, -18.0± 1.5, -4.9 ± 
0.9), (-10.1+2.1, -15.9±0.8, - 9.2±1.0),and (-9.9± 1.5, -20.9± 
0.0, -1.4 ± 0.9) from Table 1 of [Torres et al. (2008)| 
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Fig. 2. Number of y6 Pic members in a sphere of 2cr radius in the 
UVW space around each member Sigma is the velocity disper- 
sion for this association. 
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Fig. 3. Distances between TWA22AB and three moving groups 
(TWA, p Pictoris, and Tucana-Horologium moving groups) for 
the past 12 Myrs. 



Since Hipparcos proper motions were derived from obser- 
vations covering only a small time interval, we also calculated 
mean spatial motions of these associations using Tycho-2 proper 



motions ( Hoeg et al. 2000 1 and we obtained nearly the same re- 
sults: (-10.1, -17.3, -5.0), (-10.8,-16.0,-9.1), and (-9.6, -20.8, - 
0.5), respectively, for TWA, (3 Pic, and Tuc-Hor. 

To test the membership of TWA22 AB to TWA, /? Pic, and 
Tuc-Hor, based on its space motion, we applied a test with 3 
degrees of freedom using their space motion measurements. We 
find the probabilities that TWA22 AB space motion be compat- 
ible with the mean space motion of yS Pic, TWA and Tuc-Hor 
are 50%, 1% and 0.5% respectively. An alternative approach is 
to perform a k-NN analysis in the UVW space (implemented 
by I Venables & Ripley 2002, R Development Core Team^OOS i 



where we computed the distance of TWA22 to all members of 
these associations. Among the k nearest neighbors to TWA22, 
the fraction of members for a given association gives the mem- 
bership probability for that group. This k-NN analysis corrobo- 
rates that TWA22 is more likely a member of p Pic than TWA 
and Tuc-Hor Both calculations tend to reject TWA and Tuc-Hor 
as a host association for TWA22 AB. 

Our final approach is to count the number of fi Pic, TWA, 
and Tuc-Hor members (or neighbors) within a sphere of fixed 
radius in the UVW space, centered at TWA22's value and at 
each individual member of these associations. Radii are selected 
proportional to the velocity dispersions of the associations. We 
found that UVW spatial densities of TWA and Tuc-Hor mem- 
bers around TWA22 are significantly lower than the averaged 
one found in both associations. In contrary, the UVW spatial 
density of fi Pic members around TWA22 (within 2cr) is similar 
to the average density in fi Pic (see Fig. [2]). 

Membership probabilities and association member density 
around TWA22 support a possible membership in the fi Pic as- 
sociation. The membership of TWA22 in TWA cannot be firmly 
evaluated because of the paucity of high-quality kinematic data 
for most of its members. 



3.2. Trace-back 

Another way to test the membership of TWA22 AB to TWA, 
P Pic, and Tuc-Hor is to use the trace-back technique to com- 
pare the Galactic space position of TWA22 AB and TWA, p Pic, 
and Tuc-Hor backward in time. To avoid using any uncertain 
values for the solar peculiar motion (Mihalas & Binney 1981 



Dehnen & Binney 1998 , Makarov 2007) by transforming helio- 
centric velocities into LSR velocities, we decided to work in a 
reference system centered, along time, on TWA22AB instead. 
We present, in Figure |3] the distance between TWA22 AB and 
the center of each association in time. We note that TWA22 is 
always closer to p Pic than to TWA and Tuc-Hor Since the mean 
motions of the associations derived here are in good agreement 
with values from the literature, the few TWA members consid- 
ered here may not have seriously affected the result. 



4. Conclusions 

Motivated by the importance of the young, very low-mass astro- 
metric binary TWA22 AB as an important calibration point for 
stellar theoretical calculations, we measured its precise trigono- 
metric parallax (57.0 ± 0.7 mas), proper motions (-175.8 ± 0.8 
mas/yr,-21.3 ± 0.8 mas/yr), and radial velocity (Vrad - 14.8 ± 
2.1 km/s). These parameters are fundamentals for determin the 
physical properties of the tight binary system (Bonnefoy et al. 
2009). 

Our high-quality astrometric measurements along with 
HARPS radial velocity measurement allow us to discuss the 
membership of TWA22AB to nearby associations. Our kine- 
matical study shows that membership by TWA22 AB in known 
young, nearby associations can be excluded except for the p 
Pictoris and TW Hydrae associations. Membership probabilities 
based on the system space motion or the use of trace-back tech- 
nique also support possible membership of TWA22 AB in the 
P Pictoris association. Membership in the TWA cannot be fully 
excluded because of the current lack of precise parallax mea- 
surements for most of its members. Our results are, to some ex- 
tent, inconclusive about the membership of TWA22 AB in TWA 
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or yS Pic, but they are consistent with that from an age analysis. 
The location of TWA22 AB on a color-magnitude diagram sup- 
ports its age being about lOMyr but cannot be determined pre- 
cisely enough to distinguish from 8 (TWA age) and 12 Myr (fi 
Pic moving group age). Precisely known trigonometric distances 
of many more TWA members, an aim of our larger astrometric 
program of observing all known TWA members, should improve 
the situation soon. 
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